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ABSTRACT 

We report on the discovery of a periodic modulation in the bright supersoft X-ray source XMMU 
J004252. 5+411540 detected in the 2000-2004 XMM-Newton observations of M31. The source exhibits 
X-ray pulsations with a period P~217.7 s and a quasi-sinusoidal pulse shape and pulsed fraction 
~7-ll%. We did not detect statistically significant changes in the pulsation period on the time 
scale of 4 years. The X-ray spectra of XMMU J004252. 5+411540 are extremely soft and can be 
approximated with an absorbed blackbody of temperature 62-77 eV and a weak power law tail of 
photon index T ~1. 7-3.1 in the 0.2-3.0 keV energy band. The X-ray properties of the source and 
the absence of an optical/UV counterpart brighter than 19 mag suggest that it belongs to M31. 
The estimated bolometric luminosity of the source varies between ~2xl0 38 and ^8xl0 38 ergs s^ 1 at 
760 kpc, depending on the choice of spectral model. The X-ray pulsations and supersoft spectrum 
of XMMU J004252. 5+411540 imply that it is almost certainly an accreting white dwarf, steadily 
burning hydrogen-rich material on its surface. We interpret X-ray pulsations as a signature of the 
strong magnetic field of the rotating white dwarf. Assuming that the X-ray source is powered by disk 
accretion, we estimate its surface field strength to be in the range 4 x 10 5 G < B n < 8 x 10 6 G. 
XMMU J004252. 5+411540 is the second supersoft X-ray source in M31 showing coherent pulsations, 
after the transient supersoft source XMMU J004319. 4+411758 with 865.5 s pulsation period. 

Subject headings: galaxies: individual (M31) — novae, cataclysmic variables — X-rays: binaries — 
X-rays: stars 
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1. INTRODUCTION 

Luminous supersoft X-ray sources (SSSs) (Kahabka & 
van den Heuvel 2006, and references therein) were first 
discovered in the Magellanic Clouds with the Einstein 
observatory and later were established as a major new 
source class based on the results of ROSAT observations. 
SSSs have very soft spectra typically described by black- 
body models with temperatures of ~ 20 — 80 eV with no 
strong hard component, and luminosities of ~ 10 35 — 10 38 
ergs s . Although supersoft X-ray sources are not a ho- 
mogeneous class, the observed properties of the major- 
ity of SSSs are consistent with those of accreting white 
dwarfs (WD) in binary systems that are steadily or cycli- 
cally burning hydrogen-rich matter (van den Heuvel et 
al. 1992). The required accretion rates in these systems 
can be as high as 10 -7 M year . Another subclass 
of SSSs are single highly evolved stars on their way to 
WD phase. In addition, a number of more luminous 
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10 ergs s 1 ) X-ray sources also classified 



as SSS have been recently discovered in nearby galaxies, 
with some of them proposed as intermediate mass black 
hole (IMBH) candidates (Fabbiano 2006, and references 
therein) . 

The nearby giant spiral M31 presents an excellent 
opportunity to study various X-ray source populations. 
Earlier observations of M31 with the ROSAT satellite re- 
vealed a significant population of supersoft X-ray sources 
(Supper et al. 1997, 2001; Kahabka 1999; Greiner 2000). 
The advent of a new generation of X-ray telescopes 
(Chandra and XMM-Newton) has allowed us to study 
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SSSs in a much greater detail (Osborne et al. 2001; 
Trudolyubov et al. 2001, 2005; DiStefano et al. 2004; 
Greiner et al. 2004; Orio 2006). Most of the SSS detected 
in M31 appear to be transient /recurrent or highly vari- 
able in X-rays. The observations with XMM-Newton and 
Chandra are starting to provide valuable information on 
the short-term variability of M31 SSSs: timing studies of 
M31 X-ray sources have led to the discovery of 865 s pul- 
sations in the transient SSS XMMU J004319.4+411758 
(Osborne et al. 2001), and significant short-term vari- 
ations in other SSSs (Orio 2006; Trudolyubov, Pried- 
horsky & Cordova 2007). 

The X-ray source XMMU J004252. 5+411540 was dis- 
covered in the M31 field by the Einstein observatory 
(source #69 in Trinchieri & Fabbiano (1991)) and de- 
tected in subsequent observations with ROSAT (source 
#58 in Primini, Forman & Jones (1993)), Chandra 
(source r2-12 in Kong et al. (2002)) and XMM-Newton 
(source #352 in Pietsch et al. (2005)). The Chandra 
and XMM-Newton spectroscopy of the source revealed a 
thermal spectrum with effective temperature of ~ 60 — 70 
eV (DiStefano et al. 2004; Orio 2006; Trudolyubov, 
Priedhorsky & Cordova 2007), placing it in the super- 
soft source class. Assuming the distance of 760 kpc (van 
den Bergh 2000), the estimated unabsorbed luminosity 
of J004252. 5+411540 in the 0.3-1.5 keV energy band can 
be as high as ~5xl0 38 ergs s -1 , making it the brightest 
persistent supersoft source in M31. 

In this paper, we report on the discovery of the 
coherent 217.7 s pulsations in the flux of XMMU 
J004252. 5+411540, using the archival data of XMM- 
Newton observations. We also study X-ray spectral prop- 
erties of the source, search for its optical/UV counter- 
parts and discuss its nature. 
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2. OBSERVATIONS AND DATA REDUCTION 

In our analysis we used the data of four longest avail- 
able 2000-2004 XMM-Newton observations of the central 
region of M31 (Table 1) with three European Photon 
Imaging Camera (EPIC) instruments (MOS1, MOS2 and 
pn) (Turner et al. 2001; Strueder et al. 2001), and the 
Optical Monitor (OM) telescope (Mason et al. 2001). We 
reduced XMM data using XMM-Newton Science Analysis 
System (SAS v 7.0.0) 3 . We performed standard screening 
of the original X-ray data to exclude time intervals with 
high background levels, applying an upper count rate 
threshold of 20% above average background level. The 
standard SAS tool barycen was used to perform barycen- 
tric correction on the original EPIC event files used for 
timing analysis. 

We generated EPIC-pn and MOS images of the source 
field in the 0.3-7.0 keV energy band, and used the 
SAS standard maximum likelihood (ML) source detec- 
tion script edetecLchain to detect point sources. We 
used bright X-ray sources with known counterparts 
from 2MASS and USNO-B catalogs (Cutri et al. 2003; 
Monet et al. 2003) and Chandra source lists to correct 
EPIC image astrometry. The astrometric correction was 
also applied to the OM images, using cross-correlation 
with 2MASS survey catalog. After correction, we esti- 
mate residual systematic error in the source positions to 
be of the order 0.5 - 1" for both EPIC and OM. 

To generate EPIC-MOS source lightcurves and spec- 
tra, we used a circular region of 18" radius centered at 
the position of XMMU J004252.5+411540. Due to the 
source proximity to the edge of EPIC-pn CCD in obser- 
vations ## 1,2,3 the source counts were extracted from 
the elliptical regions with semi-axes of 18" and 15", in- 
cluding more than ~ 70% of the source energy flux. The 
adjacent source-free regions were used to extract back- 
ground spectra and lightcurves. The source and back- 
ground spectra were then renormalized by ratio of the 
detector areas. In this analysis we use valid pn events 
with pattern 0-4 (single and double) and pattern 0-12 
(single-quadruple) events for MOS cameras. Because of 
the extreme softness of the source X-ray spectrum, only 
data in the 0.2-3 keV energy range was used in the spec- 
tral analysis. For timing analysis, we used data in the 
0.2 — 1 keV energy band to maximize sensitivity. To syn- 
chronize both source and background lightcurves from in- 
dividual EPIC detectors, we used the identical time filter- 
ing criteria based on Mission Relative Time (MRT), fol- 
lowing the procedure described in Barnard et al. (2007). 
The background lightcurves were not subtracted from the 
source lightcurves, but were used later to estimate the 
background contribution in the calculation of the source 
pulsed fractions. 

The energy spectra were grouped to contain a mini- 
mum of 25 counts per spectral bin in order to allow \ 2 
statistics, and fit to analytic models using the XSPEC 
v.12 4 fitting package (Arnaud 1996). MOS1 and MOS2 
data were fitted simultaneously, but with normalizations 
varying independently. For timing analysis we used stan- 
dard XANADU/XRONOS v.5 5 tasks. ' 

In the following analysis we assume M31 distance of 

3 See http://xmm.vilspa.esa.es/user 

4 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 

5 http://heasarc.gsfc.nasa.gov/docs/xanadu/xronos/xronos.html 



760 kpc (van den Bergh 2000). All parameter errors 
quoted are 68% (la) confidence limits. 

3. RESULTS 

Based on the Chandra aspect solution, limited 
by ~ 0.5" systematics, the source location is 
a = 00 /l 42 m 52.52 s , 5 = 41°15'40" (2000 equinox) 
(Kong et al. 2002; Voss & Gilfanov 2007) (Fig. 1, left). 
We used the data of XMM-Newton/ 'OM observations to 
search for UV counterparts to the source. We did not 
detect any stellar-like objects in the OM images down 
to the limit of - 19 m in the OM UVW1 (291 nm) band 
within the error circle of XMMU J004252.5+411540 (Fig. 
1, right). 

3.1. X-ray Pulsations 

We performed timing analysis of XMMU 
J004252. 5+41 1540 using the data from all three 
XMM-Newton/EPIG detectors in the 0.2-1 keV energy 
band. After a barycentric correction of the photon 
arrival times in the original event lists, we performed a 
Fast Fourier Transform (FFT) analysis using standard 
XRONOS task powspec, in order to search for coherent 
periodicities. For the analysis of XMM-Newton data, 
we used combined synchronized EPIC-pn and MOS 
lightcurves with 2.6 s time bins to improve sensitivity. 
We found strong peaks in the Fourier spectra of data 
from the 2000 Jun. 25 and 2002 Jan. 6 XMM-Newton 
observations at the frequency of ^4.6xl0~ 3 Hz (Fig. 2, 
left panels). The strengths of the peaks in the individual 
Fourier spectra (Fig. 2) correspond to the period 
detection confidence of ~ 1.5 x 10~ 3 and - 8 x 10~ 7 for 
these observations (Vaughan et al. 1994). 

To estimate the pulsation periods more precisely, we 
used an epoch folding technique, assuming no period 
change during individual observations. The most likely 
values of the pulsation period (Table 2) were obtained fit- 
ting the peaks in the x 2 versus trial period distribution 
with a Gaussian. The period errors in Table 2 were com- 
puted following the procedure described in Leahy (1987). 
Then the source lightcurves were folded using the periods 
determined from epoch folding analysis. The resulting 
folded lightcurves of the source in the 0.2-1 keV energy 
band during 2000 Jun. 25 and 2002 Jan. 6 observations 
are shown in Fig. 2 (right panels). Although the power 
density spectra of the 2001 Jun 29 and 2004 Jul. 16 ob- 
servations have a much weaker (with excess power <20) 
peaks at the pulsation frequency, we still performed a 
search for the pulsations in the data of these observation 
using epoch folding technique. The epoch folding anal- 
ysis of the source lightcurves obtained during 2001 Jun 
29 and 2004 Jul. 16 observations revealed distinct peaks 
in the \ 2 versus trial period distribution at the periods 
of 217.66 s and 217.75 s (Table 2). We did not detect 
a significant change of the pulsation period between the 
four observations used in the analysis. 

The source demonstrates quasi-sinusoidal pulse profiles 
in the 0.2-1 keV energy band during all four observations. 
The pulsed fraction, defined as (I m axT m i n )/(I max +I min ), 
where I max and I m ; n represent source intensities at the 
maximum and minimum of the pulse profile excluding 
background photons, is somewhat lower during 2001 Jun. 
29 observation (6.7±1.2%) when compared to the three 
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other observations with pulsed fractions of 9.1-10.9% 
(Table 2). 

To investigate energy dependence of the source pulse 
profile, we extracted light curves in the soft (0.2-0.5 keV) 
and hard (0.5-1 keV) bands for 2002 Jan. 6 observation, 
and folded them at the corresponding best pulsation pe- 
riod (Figure 3). As can be seen from Fig. 3, there is 
a significant difference between the pulse profiles at low 
and high energies. The peak of the pulse profile in the 
0.5-1 keV band is shifted compared to that in the 0.2-0.5 
keV band. In addition, the pulsed fraction in the 0.5-1 
keV band (12.8±2.1%) appears to be higher than that at 
lower energies (7.5±1.2%). 

3.2. Variability on the Time Scale of Hours 

The source is also variable on a time scales of hours 
during all XMM-Newton observations. The characteris- 
tic timescale of the variability varies between ^6000 s 
and ^8000 s. This variability can be either aperiodic 
or quasi-periodic: the existing observations are not long 
enough to make a definitive conclusion on its nature. The 
analysis of the long, 2001 Oct. 5 Chandra/ ACIS-S ob- 
servation reveals somewhat similar variability on a time 
scale of ~14000 s. 

3.3. X-ray Spectra 

The pulse phase averaged spectra of XMMU 
J004252. 5+411540 source are extremely soft, with only a 
small fraction of photons detected above 1 keV (Fig. 4). 
Nevertheless, the source spectra can not be adequately 
fit with a single soft spectral component (i.e. blackbody 
or other thermal models): an additional hard component 
is required to approximate the excess at higher energies 
(DiStefano et al. 2004; Orio 2006). We used a sum of the 
absorbed blackbody and power law models to approxi- 
mate EPIC-pn and MOS spectra in the 0.2-3 keV energy 
range (Fig. 4). The best- fit spectral model parameters 
are shown in Table 2. The soft blackbody component has 
characteristic temperatures kT ~62-67 eV and emitting 
radii R b b ~ 13000 — 20000 km, and estimated bolometric 
luminosity L bo i ~ (5 — 8) x 10 38 ergs s _1 . The hard power 
law component is relatively weak, contributing less than 
1% of the total flux in the 0.2-1 keV energy range. The 
photon index of the power law component is poorly con- 
strained for all four observations: the overall quality of 
the fit and soft component parameters are not partic- 
ularly sensitive to the choice of the photon index. The 
equivalent absorbing column required by these model fits, 
iVn ~ (1.1 — 1.6) x 10 21 cm -2 is somewhat higher than 
the expected foreground column in M31 direction 7 x 10 20 
cm~ 2 (Dickey & Lockman 1990). 

Although our simple two-component model provides 
good overall description of the source spectra, the rela- 
tively high reduced Xr ~ 1.2- 1.7 of the fits suggest a more 
structured shape of the source spectrum and a need for 
an additional components (see also Orio 2006). The ap- 
proximation with the sum of absorbed blackbody and 
power law models leaves bump-like residuals in the 0.3- 
0.8 keV energy range, especially evident in the higher 
quality 2002 Jan. 6 data (Fig. 4). These residuals can 
indicate the presence of discrete emission/absorption fea- 
tures intrinsic to the source, or be a result of incomplete 
spectral background subtraction (the source is embedded 



in the bright unresolved emission, that have a soft ther- 
mal spectrum rich in emission lines (Shirey et al. 2001)). 

To investigate the effect of additional discrete model 
features on the continuum spectral parameters, we added 
a Gaussian line to the existing two-component model, 
and used it to approximate EPIC-pn spectrum of 2002 
Jan. 6 observation. The addition of the Gaussian line 
at ^0.53 keV greatly improved the quality of spectral fit 
(X 2 = 141.1(143 dof) with A X 2 ~ 106 for 3 additional 
parameters), and led to a significant change in the soft 
component temperature and emitting radius: kT = 76+ 3 
eV, R bb = 6998±l||i km and an absorbing column depth, 
Nh = (9±l)xl0 20 cm~ 2 . As a result, the estimated bolo- 
metric luminosity of the source dropped from ~ 7.6 x 10 38 
ergs s _1 to <~ 2.0 x 10 38 ergs s . The application of this 
model to the data of 2000 Jun. 25 observation yields 
similar results: kT = 74t 5 3 eV, R bh = 6602±i||| km, 
Nh = (8 ± 2) x 10 20 cm -2 and bolometric luminosity of 

- 1.7 x 10 38 ergs s" 1 with xl = 115 for 111 dof. We also 
tried to add a set of the discrete absorption features (ab- 
sorption lines and edges) to the original two-component 
model. This approach produced similar results: increase 
of the blackbody model temperature kT ^ 72 — 80 eV 
and decrease of the emitting radius R bb ~ 7000 — 12000 
km, and overall decrease of the bolometric luminosity to 

- (2 - 4) x 10 38 ergs s" 1 . 

This is yet another demonstration of the huge uncer- 
tainty of SSS luminosity estimates, that has poorly un- 
derstood systematic component determined by the qual- 
ity of instrument calibration, background subtraction, 
foreground absorbing column estimation and choice of 
spectral model. This uncertainty extends well beyond 
standard statistical errors of the absorbed simple black- 
body fits, usually adopted in the spectral analysis of 
SSSs, and can affect final conclusions about their proper- 
ties. In the case of XMMU J004252. 5+411540, a simple 
two-component model fits result in the estimated source 
luminosities that are supcr-Eddington assuming 1.4 M Q 
compact object accreting hydrogen-rich material (DiS- 
tefano et al. 2004; Orio 2006, this work). Recently, 
Orio (2006) performed spectral analysis of XMM/EPIC- 
MOS data on XMMU J004252. 5+41 1540 and obtained 
the lower limit on the source bolometric luminosity of 
9.5 x 10 38 ergs s" 1 for 2001 June 29 and 2002 Jan. 6 ob- 
servations, and suggested that the source could be linked 
to the ultraluminous supersoft X-ray sources found in 
other nearby galaxies and interpreted as IMBH candi- 
dates (Kong & DiStefano 2005). Our analysis of the same 
datasets yields consistently lower source luminosity esti- 
mates (~ (5 — 8) x 10 38 ergs s _1 ) for the same model 
approximation, that become even lower (~ 2 x 10 38 ergs 
s^ 1 ) when a more complex models are used for spectral 
fitting. 

4. DISCUSSION 

The X-ray pulsations and supersoft spectrum of 
XMMU J004252. 5+41 1540 imply that it is almost cer- 
tainly an accreting white dwarf in a binary system. The 
217.7 s pulsation period is the shortest known among 
SSSs to date. It is too short to be interpreted as bi- 
nary orbital period even for a double-degenerate sys- 
tem. The remaining possible explanations for the ob- 
served modulation include stellar rotation and non-radial 
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g-mode pulsations of the white dwarf. Periodic varia- 
tions on a time scales of 1000-2500 s have been recently 
detected in the postnova SSSs (Drake et al. 2003; Ness 
et al. 2003) and in one of the prototype SSS, CAL 83 
(Schmidtke & Cowley 2006), and interpreted as a sig- 
nature of non-radial pulsations of the white dwarf. In 
the case of XMMU J004252.5+411540, the short period, 
long-term stability of the modulation and the absence 
of multi-periodicity argue against that interpretation. 
Therefore, it is more plausible to assume that the ob- 
served modulation results from rotation of a magnetized 
accreting white dwarf. The strength of the WD mag- 
netic field should not be extremely high, so its spin and 
orbital periods are not locked. The energy spectrum of 
XMMU J004252. 5+41 1540 is typical for luminous super- 
soft X-ray sources and soft intermediate polars. Assum- 
ing the source is located in M31, its observed luminosity 
(~ 3 x 10 37 ergs s _1 ) is too high for intermediate polars 
(Haberl & Motch 1995). On the other hand, the absence 
of a bright optical counterpart makes the interpretation 
of the source as Galactic foreground system unlikely. The 
remaining possibility is that XMMU J004252. 5+41 1540 
is a luminous supersoft source located in M31. In that 
case, XMMU J004252.5+411540 is the second M31 SSS 
with coherent X-ray pulsations detected, after the super- 
soft transient XMMU J004319.4+411758 showing pulsa- 
tions with 865.5 s period (Osborne et al. 2001; King et 
al. 2002). 

The results of our spectral analysis suggest that the 
estimated source luminosity is probably close or may 
slightly exceed the isotropic Eddington luminosity for 
a typical white dwarfs. Assuming that the source lu- 
minosity (~ 2 x 10 ergs s _1 ) is produced in the pro- 
cess of burning of the hydrogen-rich material on the WD 
surface, we can estimate the required accretion rate as 
M ~ 5 x 10~ 7 M Q year" 1 or ~ 3.2 x 10 19 g s" 1 (assum- 
ing the energy yield of nuclear burning 0.007Mc 2 ). Such 
a high accretion rate would imply either Roche lobe ther- 
mal time scale mass transfer from a near-main-sequence 
companion with a mass larger than that of the WD, or a 
mass transfer in a wide symbiotic (WD+red giant) sys- 
tem (van den Heuvel et al. 1992; Kahabka & van den 
Heuvel 2006). 

If the X-ray pulsations observed in XMMU 
J004252. 5+411540 result from rotation of a mag- 
netized WD, accreting from the Keplcrian disk (King 
et al. 2002), one can estimate the strength Bq of the 
WD magnetic field. Assuming a dipole configuration of 
the pulsator magnetic field, the magnetosphcric radius 
inside which the accretion flow is channeled by the field 
is (Davidson & Ostriker 1973; Lamb, Pethick & Pines 
1973; Ghosh & Lamb 1979a) 
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where M denotes the accretion rate through the disk, 
Mx is a WD mass, /i = B^R^ is the magnetic moment 
of the WD, Mi 9 = M/10 19 g s" 1 , m x = M x /M , and 
M34 = m/10 34 G cm 3 . The matter will be magnetically 
channeled to the surface, if Rm > -R w d- This condition 
is satisfied for a surface field 

B > 2.4 x 10 5 i?g 5/ X /2 mj/ 4 G (2) 



where Rg = i? w d/10 9 cm. For a 1M white dwarf with 
Rwd ~ 10 9 cm accreting matter at a rate M ~ 3.2 x 10 19 
gs-\ S >4x 10 5 G. 

In order to keep the "propeller" effect from prevent- 
ing effective accretion onto white dwarf, the magneto- 
spheric radius, Rm must be less than corotation radius 
R c = (GMj/^ 2 ) 1 / 3 ? 2 / 3 , where P is the WD spin pe- 
riod (Illarionov & Sunyaev 1975). As a result, we obtain 



Bo < 8.6 x 10 3 i? 9 - 3 M 1 1 9 /2 m x /6 P 7 / 6 G 



(3) 

For the observed period P = 217.7 s, mx = 1, Rg ~ 1 
and M\g ~ 3.2, the estimated white dwarf surface mag- 
netic field Bo < 8 x 10 6 G. 

In the highly magnetized SSS scenario the accreted 
matter is channeled by the field and steadily arrives at 
the magnetic poles. Ideally, most of the hydrogen burn- 
ing is expected to occur near the poles (King et al. 2002), 
creating bright hot spots. The observed energy depen- 
dence of the pulse profile supports this picture, suggest- 
ing a higher temperature of the emitting region respon- 
sible for the modulation. The estimated fractional polar 
cap area / ~ (R wd /2R M ) ~ {R wd /2R C ) ~ 0.09 (David- 
son & Ostriker 1973) implies a relatively large size of 
the accretion (hot spot) region, in agreement with quasi- 
sinusoidal pulse profile and low amplitude of the modu- 
lation (~10%). 

In the accretion process, angular momentum is trans- 
ferred to the white dwarf and its rotation period should 
change (Pringle & Rees 1972; Rappaport & Joss 1977; 
Ghosh & Lamb 1979b). Assuming that disk accretion is 
occurring in the system, and that all the specific angular 
momentum of the accreting material is transferred to the 
white dwarf at the magnetospheric radius Rm ^$ Rc, one 
can estimate maximum rate of change of the spin period 
as 



P 7 x 10~ 8 Mi 9 /, 



50 TO X 



!/3p7/3 



s yr 



(4) 

where /50 is the moment of inertia of the white dwarf 
in units of 10 50 g cm 2 . For the typical white dwarf pa- 
rameters (-Z50 ~ 1, mx = 1), the observed spin period 
of P = 217.7 s, and M\g ~ 3.2, the estimated spin-up 
rate is \P\ <0.065 s yr -1 . The period change of such 
magnitude could, in principle, be detected in XMMU 
J004252.5+41 1540 with sufficiently long (> 50 ks) XMM- 
Newton observations spaced by a period of > 5 years 
(Table 2). The measurement of the period change can be 
than used to further constrain the accretion rate in the 
system, providing an independent, better estimate of the 
bolometric luminosity of the source than presently pos- 
sible with spectral analysis. Unfortunately, only one of 
the existing observations (2002 Jan. 6) is sufficiently long 
to obtain good estimate of the pulsation period. Future 
long XMM-Newton observations of the central region of 
M31 would allow a more detailed study of the evolution of 
the pulsation period of XMMU J004252. 5+411540, and 
possibly improve our understanding of this and other lu- 
minous supersoft X-ray sources, both Galactic and ex- 
tragalactic. 
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TABLE 1 



XMM-Newton Observations of M31 Used in the Analysis. 


Obs.# 


Date, UT Date, TJD 


Obs. ID Mission/Instrument Mode/ 


RA (J2000) a Dec (J2000) a Exp.(MOS/pn) b 








Filter 


(h:m:s) 


(d:m:s) 


(ks) 


1 


2000 Jun. 


25 


0112570401 XMM-Newton/EPIC Full/Medium 


00:42:43.00 


41:15:46.1 


28.9/24.9 


2 


2001 Jun. 


29 


0109270101 XMM-Newton/EPIC Full/Medium 


00:42:43.00 


41:15:46.1 


29.0/24.9 


3 


2002 Jan. 


06 


0112570101 XMM-Newton/EPIC Full/Thin 


00:42:43.00 


41:15:46.1 


63.0/49.9 


4 


2004 Jul. 


16 


0202230201 XMM-Newton/EPIC Full/Medium 


00:42:42.12 


41:16:57.1 


19.6/16.4 



a pointing coordinates 

b instrument exposure used in the analysis 
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TABLE 2 

X-ray Pulsation Parameters and Blackbody Spectral Fit Information for XMMU J004252. 5+411540. 



Timing Parameters (BBODYRAD+POWERLAW)*WABS Spectral Model Parameters 

Obs. Period PF . 2 -i ke v N H kT bh R hh Photon Flux b L x c ix d X 2 Instrument 





(s) 


(%) a 


(xlO 20 cm" 2 ) 


(eV) 


(km) 


Index 








(d.o.f) 




1 


217.60(23) 


10.9 ±2.0 


12 ± 1 


65± 2 


149031™ 


2-59li;- 


4.66 ±0.05 


3.22 


3.1 


135.8(114) 


pn 








11 ± 1 


67 ±2 


12912ll?S 


1 QO+1-28 


5.02 ±0.06 


3.47 


2.7 


87.2(79) 


MOS 


2 


217.66(32) 


6.7 ± 1.2 


14 ±1 


63 ±2 


18240t 2 2 5 7 g 


2-63±i"i| 


4.17 ±0.06 


2.88 


3.8 


133.1(114) 


f 








12 ± 1 


63 ±2 


1 '7'70'7"l~2566 
1 ' ' °' -2664 




5.01 ±0.08 


3.46 


3.7 


85.8(76) 


MOS 


3 


217.73(08) 


9.1 ± 1.2 


16 ±1 


66 ± 1 


19553™ 


i yir+0.67 


5.17 ±0.05 


3.57 


5.5 


247.3(146) 


pn 








13 ±1 


63 ±1 


18462±^ 


9 QR+0.64 
^- OU -0.69 


4.78 ±0.06 


3.31 


3.8 


137.2(96) 


MOS 


4 


217.75(48) 


10.7 ±2.5 


12 ± 1 


62± 3 2 


17569tg 2 ? 


3.06±i;?| 


4.58 ±0.08 


3.16 


3.4 


134.5(100) 


pn 



a pulsed fraction in the 0.2 — 1 keV energy band, defined as (Imax-Imin)/(Imax+I m in): where l m ax and I m i n represent source background- 
corrected intensities at the maximum and minimum of the pulse profile 
b absorbed model flux in the 0.2 — 1 keV energy range in units of 10 — 13 erg s _1 cm~ 2 

c absorbed luminosity in the 0.2 — 1 keV energy range in units of 10 37 erg s _1 , assuming the distance of 760 kpc 

d cstimated unabsorbed luminosity in the 0.2 — 1 keV energy range in units of 10 38 erg s -1 , assuming the distance of 760 kpc 

f part of the source flux falls into the gap between EPIC-pn CCDs 



M31 Center/XMM-Newton/EPIC-MOS (0.3-7 keV) M31 Center/XMM-Newton/OM (UVW1 filter) 




30.0 0:43:00.0 30.0 42:00.0 55.2 52.8 50.4 0:42:48.0 

R. A. (J2000) R. A. (J2000) 



Fig. 1. — Left: Combined 0.3-7 keV XMM/EPIC-MOS image covering central region of M31. The position of supcrsoft pulsating source 
XMMU J004252. 5+411540 is marked with an arrow. Right: XMM-Newton/OM UVW1 (291 nm) band image of M31 field taken on Jan. 
6, 2002. The image is a 2' X 2' square centered on the XMMU J004252.5+411540 position. The localization of XMMU J004252.5+411540 
is shown with black circle of 1.5" radius (3<r). 
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POWER DENSITY SPECTRUM 
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Fig. 2.— (Left) Power spectra of XMMU J004252.5+411540 obtained using the data of 2000 Jun. 25 (upper panel) and 2002 Jan. 6 
XMM-Newton/EPIC observations (lower panel) in the 0.2-1 keV energy band. (Right) Corresponding pulse profiles folded with most likely 
pulsation periods, corrected for background. 
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Fig. 3. — Background-corrected normalized X-ray lightcurves of XMMU J004252. 5+411540 during the 2002 Jan. 6 observation folded 
at the corresponding best period (Table 2) in the 0.2-0.5 and 0.5-1 keV energy bands (upper and middle panels) along with their hardness 
ratio (bottom panel). 
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Fig. 4. — EPIC-pn count spectra and model ratios of the source during during four XMM-Newton observations used in the analysis. The 
best-fit absorbed blackbody plus power law model approximation (Table 2) is shown with thick histograms. 



